Application of a renormalization group algorithm to nonequilibrium cellular automata 

with one absorbing state 
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We improve a recently proposed dynamically driven renormalization group algorithm for cellular 
automata systems with one absorbing state, introducing spatial correlations in the expression for the 
transition probabilities. We implement the renormalization group scheme considering three different 
approximations which take into account correlations in the stationary probability distribution. The 
improved scheme is applied to a probabilistic cellular automaton already introduced in the literature. 



I. INTRODUCTION 

In order to study cellular automata systems displaying 
a second order irreversible phase transition to an absorb- 
ing state characterized by a scalar order parameter, a 
dynamical renormalization group (RG) algorithm ffl has 
been recently proposed. These cellular automata models 
are in the directed percolation (DP) universality class. 
The method is based on a Dynamically Driven Renor- 
malization Group (DDRG) scheme (for a recent review 
about DDRG see Q ) , which has been successfully applied 
to self-organized critical phenomena, as sandpile models 
H and forest fire models Q. 

The basic idea introduced in [Q is to couple a real 
space RG scheme to a stationary condition that drives 
the RG group equations through the parameter space. 
The stationary equations, involving the stationary dis- 
tribution, have to be approximated, since the form of the 
stationary probability distribution is not known a priori, 
as in the case of systems in equilibrium. Oliveira and Sat- 
ulovsky |ffl] showed, as proposed in Q , that results can be 
improved using more refined approximations for the sta- 
tionary probability distribution. The expression for the 
transition probability used in ]lj consists in a product 
of independent one-site transition probabilities at every 
step of the RG transformation. 

In this work we exploit another aspect of the scheme, 
not considered before, in order to include additional cor- 
relations. In fact, correlations can be also introduced 
in the renormalization scheme if we allow the transition 
probability to depend upon more neighbors. In the case 
of non-equilibrium models, the space in which the RG 
flows is the space spanned by the transition probabili- 
ties. As we will see later, our approach broadens this 
space providing more degrees of freedom to the RG tra- 
jectories that flow towards the fixed points. 

We apply the modified RG scheme to a probabilistic 
cellular automaton (PCA) with one absorbing state al- 
ready introduced in the literature M. Using a block 



renormalization that properly treats the nature of the 
absorbing state, we figure the value of the critical expo- 
nent of the divergence of the spatial correlation length, 
v± , using three different approximations involving corre- 
lations among clusters of 1, 2, and 4 neighboring sites in 
the lattice. Our calculations for i>±_ , at small orders of ap- 
proximation in the mean field scheme for the stationary 
distribution, give better values than the ones reported in 
reference [[lj. 

The paper is organized as follows. We begin with a 
brief description of the model proposed in |5| and the gen- 
eral renormalization scheme. After this, we define the al- 
gorithm used in this work and present the values obtained 
for v± . Finally, after mentioning the ideas involved in the 
simulation technique used to study the model, we present 
the values of the whole set of critical exponents for the 
PCA obtained by means of dynamical numerical simula- 
tions and stationary simulations. Our results are in well 
agreement with the values corresponding to (1 + 1) DP, 
and differ considerably from the ones reported in refer- 
ence pj. 



II. THE MODEL 

The model studied in j^] is a one-dimensional cellu- 
lar automaton in which each site can be either vacant, 
<7j = 0, or occupied by a particle, <7j = 1. At each 
time step, the state of a given site will depend only on 
its previous state and the previous state of its nearest 
neighbors. The transition probability T(a\a J ) from state 
a' = (a[, <j' 2 , • ■ • , a' L ) to state a — (at, a%, ■ ■ ■ , a£) will be 
given by the product 



T(a\a') = J\_T(ai\a' i _ x ,a' i ,a' i+x ), 



(1) 



where L is the number of sites and r(ai\a' i _ 1 , a[, a' i+1 ) is 
the one-site transition probability given by the following 
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(2) 



This probabilistic cellular automaton models a transi- 
tion from elementary rule 4 to elementary rule 22 (follow- 
ing Wolfram's nomenclature scheme) Q . The system has 
two critical points, one of them at p± — and the other at 
Pi ~ 0.75. The critical point at p\ — will be shown to 
be trivial, in contrast to the result obtained in reference 
||. Consequently, we will apply the RG to study the 
non-trivial transition point. Rule r(0|000) = 1 implies 
that the vacuum state is indeed an absorbing state. 



III. RENORMALIZATION SCHEME 

The RG scheme proposed in Q is a real space RG 
scheme jllfl in which one renormalizes the transition 
probability, T. The RG flow takes place in the space of 
parameters defining T. A blocking procedure transforms 
cells of b sites into one site at the new scale. In order 
to account for the fact that the vacuum state is absorb- 
ing, a cell devoidcd of particles will always renormalizc 
into an empty site. Cells with at least one particle have 
been chosen to renormalize into an occupied site. Other 
options have been tried, but they do not preserve the 
existence of the absorbing state. 

Let a = (oi, (72, • • • , ox) be the state of a system with L 
degrees of freedom, and the vector S = (Si,S2, ■ • • , Sl<) 
be the state of the renormalized system with V = L/b 
degrees of freedom, where b is the size of the renormal- 
ization block. 

The conditional probability of state S given state a, 
lZ(S\a), must satisfy 



TZ(S\a) > 0, 



(3) 



Given T and the probability of a state a' at time t, 
W(cr'), one can write the joint probability of state a' at 
time t and state a at n time steps later, W n (a, a'), by 
simply applying T to W(a') n successive times 



W n (a,a')=T n (a\a')W(a'). 



(4) 



In addition, in the stationary regime, the probability dis- 
tribution W(a) must satisfy 



W{o)=Y^T n {a\a')W{a'), 



(5) 



for any value of n. 

In the same way one can write these expressions at the 
coarse grained level. Denoting by T(S, S') the probabil- 
ity of occurrence of state S' at a given time and state S 
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m 



one time step later, the RG transformation is obtained 
imposing 



w(s, s') = j2Y, nsw)ns'W)w n {o, a'), (&) 



a a' 

from which follows 



W(S') = W(S, S') = J2 n(S'\a')W(a>). (7) 

S a' 

Once knowing the transition probabilities at the coarse 
grained scale, one can easily build the rescaled transition 
matrix, T(S\S'), as 



T(S\S') = 



W(S, S') 
W(S') 



(8) 



Using now equations (Q), (^|), and (0), we obtain the 
final expression for the renormalization equations 

rp, q \ q n H^nSW)n{S'\a')T^a\g')W{a') 

T{blS) ~ E„, rc(sV)wV) ■ () 

Let us note that, while this equation expresses rescaled 
transition probabilities, T, in terms of transition proba- 
bilities at a lower scale, T, the stationary weight of each 
state present in equation (j^), W(a'), is yet unknown. 
Contrary to the case of closed systems in thermal equi- 
librium, we do not know a priori the expression for the 
stationary probability distribution. However, including 
equation (||) , one can get a closed set of equations to solve 
at each renormalization step. The stationarity condition 
(H) is actually essential in driving the RG equations (||) 
through parameter space. Equation (||), together with a 
given approximation for the stationary probability, pro- 
vide then a well defined RG transformation T — > T. 

In practice, equation (||) can hardly be solved and one 
must resort to approximations. The values of critical 
exponents obtained using the present RG approach are 
expected to improve as these approximations improve. 
We have used three different levels of approximations, in 
which correlations among clusters of 1, 2, and 4 sites are 
considered respectively. 

In this work we will be concerned with another way to 
improve the RG scheme. An important point is that, in 
order to solve equations (|^), an assumption needs to be 
made on how the transition probability between states 5" 
and S, T(S\S'), depends on local transition probabilities 
at the coarse grained level. This choice will determine 
the degree of proliferation that the RG will have, since 
the form of the renormalized transition probability will 
be preserved along the RG trajectories. 

In the former approach [jj], the authors carried out 
the most general RG transformation considering one-site 
transition probabilities. By preserving the form of the 
renormalized transition probabilities, the RG trajecto- 
ries are found to flow to the attractive fixed points in a 
five-dimensional space spanned by the dynamical param- 
eters. In contrast to usual RG methods in which new 



couplings arise at each step of the transformation, this 
RG procedure is not able to proliferate the dynamical pa- 
rameters since the form of the transition probabilities are 
kept fixed at the coarse grained level. However, new dy- 
namical parameters can be considered from the very be- 
ginning if we allow the transition probabilities to depend 
upon more neighbors, i.e., introducing more correlations 
in them. This feature of the method should be compared 
with other dynamical RG procedures [Q where the in- 
troduction of new couplings since the very beginning is 
an alternative way to carry out the RG transformation. 

We propose, then, a form for the coarse grained tran- 
sition probability consisting of a product of independent 
two-site transition probabilities instead of one-site tran- 
sition probabilities. Denoting the position of each lattice 
site with i, the new transition probability is defined at 
even time steps as 

L/2 

T(S\S')= Yl ^(^Si+ilSl^^^S^S'^), ( 10 ) 

k=l,{i=2k) 

and at odd time steps as 

L/2 

T{S\S')= JJ r(Si, Si + i | S'i, S' i+1 , S' i+2 ). (11) 

fc=l,(i=2fc+l) 

Here we have used the same symbol, r, to indicate a 
different type of transition probability than the ones ap- 
pearing m ([!]). In formulas @ and (y) 

periodic bound- 
ary conditions are assumed. 

One can retrieve at any time one-site transition prob- 
abilities knowing both two-site probabilities and the sta- 
tionary distribution. It is straightforward to show that 

r{Si\S i _ l , S i+1 ) = y \ (Sj, Si^i\S i _ 1 , Si, S^+i, $1+2) 

Si+i,S' i+2 

W(S' i _ 1 , Sj, S i+i , S i+2 ) 



R{SkW2k-l,CT2k) > 0, 



W(S i _ 1 ,S' i , S' i+ll 



(12) 



Using expressions ([K]) and (y) in equation (g) one 
can implement the RG transformation. The algorithm 
we used is explained in more detail in the next section. 



IV. RENORMALIZATION ALGORITHM 

We have used a temporal coarse graining of two time 
steps (n — 2). The blocking operator, TZ, was chosen in 
the same way as in (jjj, renormalizing cells of size b = 2 
into one site 



L/2 



(13) 



and 



(14) 



(15) 



In order to preserve the nature of the absorbing state, we 
have also required R to satisfy 



and 



fl(G|0,0) = 1, 



R{0\a 2k -i,(r2k) = 0. 



(16) 



(17) 



whenever a 2 k-i ^ 0, or <j 2 k ^ 

The diagram in figure (|l|) indicates how two-site tran- 
sition probabilities are renormalized. Indexes appearing 
in equation (18)- (EM refer to this diagram. 



k=l 



8, S 2 
FIG. 1. Diagram showing the blocking scheme proce- 
dure. Numbers correspond to the indexes used in equa- 
tion 

Using equations (||), (|i"o|), and (pd|), we can write down 
the expression relating r to r, which is given by 

t(S%, S 2 , 15*3, S4, S5, Sq) = [N(Sa, S4, S5, Sq)] 

E R(Si\ax, cj 2 )R{S 2 \o 3 , 0-4)^(53^11, CT12) 

(TiCT2(T3Cr4Crii---CTi8 

R(S 4 \cri3, <Ji4,)R(S 5 \ai 5 ,a 16 )R(S 6 \an, a ls ) 
D(ai,a 2 ,a 3 , os\a Y \, • • • , <J\%)W {<J\\, ■ ■ ■ , <7 l8 ), (18) 

where 

D{a\a 2 , 03, <T4|(Tii, (T12, CT13, CT14, fJi 5 , CT16, CT17, cr l8 ) = 

E ^( CTl ' ^alfg, 06, 0"7, Cr 8 )f(f7 3 , (T4 I (Tj , (Tg , (7g , (Tio) 

tT5'"(T10 

t((T 5 , CT 6 |(Tii, CTi2, f7i3, (Ti4)r(cr7, <T8\^13, ^15, &w) 

t(o"9, Crio|cri5, CTi 6 , (T17, (Tis) , (19) 



with 



and 



JV(S 3 , S4, ^5, Sq) = ^ R(Ss\(Tn, o*i2)-R(iS4|o"i3j ffu) 

(Til -"CIS 

i2(5 B ki 6 , ffi 8 )ii(5a|(7i7, o-ibJW^h, • ■ • , a 18 ). (20) 

Since we do not know a priori the stationary weights, 
W(au, • • • , els), we need an approximate method to es- 
timate them. The simplest approximation, sometimes 
known as simple mean field approximation, consists in 
neglecting correlations among different sites. That is 

18 

W(a 11 ,---,a 18 )=l[W(a l ), (21) 

»=n 

where W((7i) is the solution of 

W(ai) = ^ T(<7ia 2 \<T3CT4:<J5V6)W(a2) 

fT2<T3<T4(T5Cr6 

W(a 3 )W{<j 4 )W(a 5 )W{a 6 ). (22) 

Correlations, however, are actually taken into account 
in the geometrical aspects of the blocking procedure, 
leading to non classical exponents. Better approxima- 
tions can be also implemented (as a reference see [p"3|). 
We used one, two, and four-site approximations. 

Equations (|l8|)-(p0|) involve each transition probability. 
Being so many terms, we are prevented from an analytical 
determination of the fixed points of the transformation. 
So, we performed our search numerically, using initial 
values for the transition probabilities that correspond to 
the original model (table (||)). 

As a technical remark, let us say that in each iteration 
of the RG, and given a set of parameters, {r}, we need to 
solve equation (^2|) (or its analogue for two and four-site 
approximations) before the next RG step. We have done 
this by iterating the equation until reaching convergence. 



0.7484 
0.7474 

Y 

0.7464 
0.7454 




0.2325 0.2332 0.2339 0.2346 0.2353 

X 

FIG. 2. Two dimensional projection of the RG flow along 
x = ?(00|0010) and y = 7(00|1000) for the simple mean field 
approximation. 

We have found one relevant parameter. Since we are 
only dealing with stationary properties of the model, it is 
reasonable to assume that this parameter is associated to 
the divergence of the spatial correlation length and not 
to the temporal correlation length . In order to calcu- 
late the eigenvalue, A, associated to that parameter, we 
have to find the linear region of the RG transformation. 

Let us take a trajectory passing close enough to the 
unstable point, and construct a sequence of numbers 
consisting in the distance between two successive points 
along that trajectory. Now, let us call r to the ratio be- 
tween two consecutive numbers in that sequence. In the 
portion of the trajectory corresponding to the linear re- 
gion of the RG transformation (around the unstable fixed 
point), one expects to see two plateaus in the values of 
r. 

3.00 I . 1 . 1 . 1 1 



V. RG RESULTS 
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The behavior of the RG equations can be described 
as follows. For values of p which are small enough, the 
set of transition probabilities flows towards an attractive 
fixed point characterized by a lattice devoid of particles. 
Increasing p above a critical value, p cr , the flow is driven 
to another attractive fixed point, consisting of a lattice 
full of particles. The value of p cr , for each level of ap- 
proximation used, can be found in table (23). Starting 
around the critical values, the representative point of the 
parameter set spends a long time near an unstable fixed 
point before leaving towards one of the two attractive 
fixed points. A projection of the RG flow along two spe- 
cific transition probabilities is shown in figure (0). 
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FIG. 3. The ratio, r, of successive distances among two 
consecutive points for one of the trajectories shown in figure 
(g). The x coordinate denotes the precedence of each point 
along the curve. 



The first plateau corresponds to a trivial parameter, 
while the value of r at the second plateau corresponds 
to the eigenvalue A of the RG transformation. Figure 
(|J) shows an example of one of such curves for the sim- 
ple mean field approximation. So, figuring the eigen- 
value A associated to the relevant parameter, we get 
u± = In 2/ In A. The value measured numerically for the 
plot in figure @ (simple mean field approximation) is 
vx = Q.965 ± 0.001. 

To the best of our knowledge, the best value of u± is 
v±_ = 1.0972 ± 0.0005 Q], and the most accurate value 
for the critical point P2, obtained by dynamical numer- 
ical simulations (see following sections), is P2 ~ 0.7513. 
Although our value of the critical exponent v± is still 
inaccurate, it is worth at this point comparing it with 
the one obtained using the same approximation (simple 
mean-field) in In that work the authors obtained 
Ux = 0.931 ± 0.005, so that our result is a better approx- 
imation for the actual vx- Our scheme is able to take 
into account correlations in a more accurate way than 
the original one. 

By increasing the order of the approximation, results 
improve, as shown in table (p3|). The value we found us- 
ing the two-site approximation was vx — 1-013 ± 0.001, 
while the value found in the four-site approximation was 
vx = 1.015 ± 0.001, which is closer to the actual one. 
Below we show the critical value, p cri for the three ap- 
proximations, as well as the corresponding value of A and 
vx 



appr. 


Per 


A 


vx 


1 


0.639825 


2.050 


0.963 ±0.001 


2 


0.681490 


1.982 


1.013 ±0.001 


4 


0.695017 


1.979 


1.015 ±0.001 



(23) 



It should be noted that the present RG scheme leads 
to fairly good values for the exponent vx already within 
lower order mean-field approximations. This fact indi- 
cates that the introduction of new correlations in the 
transition probabilities plays a relevant role within the 
lower order mean-field approaches. For mean-field ap- 
proximations of order higher than two, the convergence 
of the scheme becomes slower. 

Reconsidering the ideas which led us to Eqs. ( [Toj ) and 
( pi] ) for the transition probabilities, one may think of 
obtaining better approximations for Vx by allowing the 
transition probabilities to depend upon even more neigh- 
boring lattice sites. While this idea is clearly right, one 
can presently not overcome, in practice, the huge amount 
of computer time needed to obtain values that are accu- 
rate enough. 



VI. SIMULATION TECHNIQUE 



A. Spreading analysis 

The basic rules governing the dynamical evolution of 
the system have been formulated in section II. As in or- 
dinary cellular automata, all lattice sites are updated si- 
multaneously. Simulations were performed on lattices of 
size L = 10000, taking periodic boundary conditions. We 
briefly discuss here the scaling theory for directed perco- 
lation which supports the spreading analysis. A detailed 
treatment can be found elsewhere jl(J. 

It should be stressed that for L finite, the steady state 
of the system is meta-stable since, due to fluctuations 
of the stochastic process, there is always a finite prob- 
ability for the system to become empty. This probabil- 
ity increases when approaching the critical point. Con- 
sequently, it is very difficult to calculate critical points 
and critical exponents by means of numerical simulations. 
Furthermore, since the transition between the stationary 
regime and the absorbing state is second order, a mean 
field treatment is not adequate. These shortcomings can 
be avoided by evaluating critical exponents related to the 
dynamic critical behavior of the system. For this purpose 
one starts, at t = 0, with a particle at the center of the 
lattice otherwise empty, i.e. a configuration very close to 
the absorbing state. Then, the following quantities are 
computed: (i) the survival probability, P(t), that is, the 
probability that at least a particle is still in the system 
at time f; (ii) the average number of particles, N(t), and; 
(iii) the average mean distance, R(t), over which particles 
have spread. Averages are taken over 5 x 10 samples, 
and runs are performed up to t — 10 4 . Finite size effects 
are avoided, since the epidemic disk never reaches the 
edge of the lattice during the simulation. Close to the 
critical point, and for long enough times, the following 
scaling laws should hold ]10| |: 

P(t) oc^^A^ii}, (24) 
N(t) ocfV{Ai 1/iy i'}, (25) 

R(t) oc t z/2 E{A i 1 /"!! }, (26) 

where A = \p — p c \, £ t = A _l/ n gives the temporal cor- 
relation length close to p c , Vu is the correlation length 
exponent (time direction), $, (f>, and S are suitable scal- 
ing functions, and 8, rj, and z are critical exponents. In 
the absorbing state, P(t) and N(t) are expected to decay 
exponentially, since correlations are short-ranged. This 
can only happen if <p(A,t) cx {A t x l v w }-^ v \\ exp(-A l/ "t) 
for t — > co. Therefore, one has from Eq. (4): 

N(t) oc A" 7 *"" exp(-A^), t -> co. (27) 

At criticality, one expects that log-log plots of P(t), 
N(t), and R(t) would give straight lines, while upward 
and downward deviations would occur even slightly off- 
criticality. This behavior would allow a precise determi- 
nation of the critical point and the critical exponents 8, 



7], and z. It should be noted that by means of Eq. ( pjj ) 
it would be also possible to calculate v\\ . 

B. Finite size scaling analysis 

As in standard second order phase transitions it is 
assumed that in the supercritical region, and close to 
the critical point, the system displays spatial correlations 
characterized by a typical length scale, £ s , which diverges 
at criticality according to 

6ocA-^, A^O, (28) 

where v± is the correlation length exponent in the spatial 
direction. The natural order parameter of the model is 
the density of particles, p, which at criticality depends 
on the system size L and A as 

p{p,L) = L-^f(AL 1 /^), (29) 

where / is a suitable scaling function and (3 is the order 
parameter critical exponent. For small positive A, and 
L —> oo, f(x) should have the following form 

f(x) cx x^ (30) 

in order to recover the well known critical behavior of the 
order parameter in the thermodynamic limit 

pocA' 3 . (31) 



VII. SIMULATION RESULTS 

Before presenting the simulation results we will briefly 
discuss the critical point at p = 0. If we start at t = 
with a random initial configuration of density po = 0.5, 
the stationary density of the system is p r* 1/8. It should 
be noted that at p = the system reaches a static sta- 
tionary state in one time step (see the evolution rules in 
table (Q)). Taking into account this observation, the sta- 
tionary density of the system can be obtained as follows. 
It is clear from the evolution rules that the probability of 
having an occupied site at t = 1 is equal to the probabil- 
ity of finding an occupied site surrounded by empty sites 
at t = 0. Since there are no correlations in the initial 
state, the stationary density of the system at p — can 
be determined as 

P = p (l-po) 2 (32) 

Then, since we used po = 0.5, the expected value of the 
stationary density is p = 1/8. 

For arbitrary small values of p, the stationary state is 
the vacuum state ||. We now consider the relaxation 
process to the vacuum state. We take as starting con- 
figuration (t — 0) any of the static stationary states at 



p = 0. Then, we follow the evolution of the density for 
p close to zero. It is clear from the evolution rules that 
after a particle is created, two particles are removed from 
the system at the next time step. Then, the system stays 
in another static state until the next creation process oc- 
curs. Suppose that a creation process happens at position 
Xi . If no new particle is created in the neighborhood of Xi 
at the next time step, correlations can not be generated 
(see the evolution rules in table @). The probability 
of two consecutive creation processes is p 2 < p. Then, 
the system can not develop long range correlations and 
a mean field analysis should be appropriate. For p close 
to zero it is possible to consider the creation process as 
a creation-induced annihilation process. We then have 

dp/dt = —pp. (33) 

Consequently, p = is a trivial critical point since the 
relaxation time behaves as r = (1/p) . This result is 
in disagreement with the one reported in reference 
probably due to a poor statistics of the simulation data. 

In the following, the results of the epidemics analysis 
at p — p2 are presented. We measure the time evolution 
of P(i), N(t), and R(t) for different values of the param- 
eter p. Log-log plots of these quantities as a function of 
time are straight lines at the phase transition and show 
curvature away from the transition. It is important to 
mention that the epidemics analysis is a very sensitive 
method since it is possible to distinguish among super- 
critical and sub-critical behavior for p values that differ 
in the fourth decimal. Our best estimation of the critical 
point is P2 = 0.7513 ± 0.0002, and the dynamical critical 
exponents are S = 0.162 ± 0.0004, rj = 0.304 ± 0.0005, 
and z/2 = 0.643 ± 0.0007. It should be remarked that 
the error bars merely indicate the statistical error ob- 
tained from regressions. The values of the dynamical 
exponents are in well agreement with those correspond- 
ing to directed percolation in (1 + 1) dimensions, as it 
was expected. 

It is possible to calculate the exponent v\\ from the 
analysis of the sub-critical behavior (see eq. (p7j)). In 
fact, the decay constant, A = ^ 1 , governing the long 
time behavior of N(t), behaves according to 

\ = fr 1 = A v *, (34) 

so, if p2 is known, we can calculate A for different val- 
ues of A. Then, a log-log plot of A vs A allows us to 
evaluate the exponent v\\ . This analysis gives an expo- 
nent v\\ — 1.738 ± 0.002 which is quite close to the value 
V\\ = 1.73 corresponding to (1+1) DP 0]. It should be 
pointed out that our value of v\\ sharply differs from the 
one reported in reference j^] {v\\ « 1.087). The error in 
the last value of the exponent v\\ is due to the fact that it 
was calculated taking into account not only sub-critical, 
but also supercritical curves. 

We have also calculated the order parameter critical 
exponent measuring the density, p, as a function of A 
in the supercritical regime (see eq.([3l|)). We obtain 



(3 = 0.277 ± 0.002, which is once again very close to 
f3 = 199/720, corresponding to (1 + 1) DP @. It should 
be mentioned that the reported value of the order param- 
eter critical exponent in reference || is (3 ~ 0.32 which 
differs around 16% from the theoretical value. This dif- 
ference may be due again to the poor statistics of the 
data. 

We finally present the finite size scaling analysis. Fig- 
ure ([|) shows a log-log plot of pL^/ u± vs AX 1 /"- 1 - for 
different values of p, and lattice sizes, X, where we have 
used (3 = 199/720 and v L = 1.0972 @ corresponding to 
(1 + 1) DP. We obtain an excellent collapse of the data 
on an universal curve, as it is predicted by eq. (p9|). 

Right before submitting this manuscript for publica- 
tion, the author of reference H improved some of his pre- 



vious results |14|. Although no new results are reported 
for the trivial critical point p\, the conclusions concern- 
ing the universality class of the model are in complete 
agreement with the ones found in this work. 



in the transition probabilities is the relevant reason for 
the improvement of the results. 

The critical exponents v±, specially for low order ap- 
proximations, are better than the ones obtained with 
schemes that make use of an independent product of one- 
site transition probabilities . 

Using very simple arguments, we have shown that 
Pi = is a trivial critical point, since the time relaxation 
constant behaves as r = (1/p) 1 . This behavior differs 
from the one reported in reference || (r = (1/p) 0,86 ). 

We have also obtained, by means of numerical sim- 
ulations, the critical point p — p2, and the whole set 
of critical exponents. The value of the critical point 
p 2 = 0.7513 ± 0.0002 is in agreement with [§. How- 
ever, we found very different values for exponents vn and 
f3. Our values of the exponents at p = pi are in well 
agreement with those corresponding to (1 + 1) DP as it 
was expected, since there is only one absorbing state for 
the system 
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VIII. CONCLUSIONS 



We have introduced a new renormalization group al- 
gorithm for probabilistic cellular automata with one ab- 
sorbing state. The new scheme introduces correlations in 
the RG procedure by allowing the transition probabilities 
to depend upon two neighboring lattice sites. Three dif- 
ferent approximations for the stationary probability dis- 
tribution have been used, namely the simple mean-field 
approximation, the pair mean-field approximation, and 
the four-site mean-field approximation. 

The new RG scheme leads to fairly good values for v± 
even within mean-field approximations of low order. This 
result shows that the introduction of spatial correlations 
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